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ABSTRACT 

Recent deep blank field submillimeter (submm) surveys have revealed a population of luminous high 
redshift galaxies that emit most of their energy in the submm. The results suggest that much of the star 
formation at high redshift may be hidden to optical observations. In this paper we present wide-area 
850-/im SCUBA data on the Hawaii Survey Fields SSA13, SSA17, and SSA22. Combining these new 
data with our previous deep field data, we establish the 850-/im galaxy counts from 2 mJy to 10 mJy 
with a > 3cr detection limit. The area coverage is 104 arcmin^ to 8 mJy and 7.7 arcmin^ to 2.3 mJy. 
The differential 850-/xm counts are well described by the function n{S) — No/{a + S^-^), where S is the 
flux in mJy, Nq = 3.0 x 10"* per square degree per mJy, and a = 0.4 — 1.0 is chosen to match the 850-/zm 
extragalactic background light. Between 20 to 30 per cent of the 850-^m background resides in sources 
brighter than 2 mJy. Using an empirical fit to our > 2 mJy data constrained by the EBL at lower fluxes, 
we argue that the bulk of the 850-^m extragalactic background light resides in sources with fluxes near 
1 mJy. The submm sources are plausible progenitors of the present-day spheroidal population. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation 



1. INTRODUCTION 

The past few years have seen a dramatic evolution in 
our knowledge of the history of star formation in optically- 
selected galaxies covering the redshift range from z = to 
z ~ 6 (see, e.g., Madau, Pozzetti, & Dickinson 1998 and 
references therein; Dey et al. 1998; Hu, Cowie, & McMa- 
hon 1998; Weymann et al. 1998; Cowie, Songaila, & Barger 
1999). However, it has long been suspected that much 
of the early star formation in galaxies takes place inside 
shrouds of dust (e.g. Partridge & Peebles 1967). While it 
is possible to attempt to correct the optical star formation 
history for the effects of dust (Heckman et al. 1998; Pettini 
et al. 1998), the corrections are far from certain. However, 
the importance of absorption and reradiation of light by 
dust in the history of galaxy formation and evolution is 
evident: the submm extragalactic background light (EBL) 
(Puget et al. 1996; Guiderdoni et al. 1997; Schlegel et al. 
1998; Fixsen et al. 1998; Hauser et al. 1998) has approx- 
imately the same integrated energy density as the optical 
EBL. Thus, an accurate reconstruction of the global star 
formation history versus redshift requires an understand- 
ing of both the rest-frame ultraviolet and the rest-frame 
far-infrared properties of the galaxy populations. 

Reradiation of stellar light by dust into the far-infrared 
produces a thermal emission peak at A 60 — 100-/im 
that is redshifted into the submm for galaxies at z > 1. 
Submm observations are unique in that the strong negative 
K-corrections at 850- /im nearly compensate for cosmologi- 
cal dimming beyond z ~ 1, thereby making dusty galaxies 
almost as easy to detect at z ~ 10 as at z ~ 1 (Blain 
& Longair 1993; McMahon et al. 1994). The revolution- 
ary new camera, SCUBA (Submillimeter Common User 
Bolometer Array; Holland et al. 1998), on the 15-m James 
Clerk Maxwell Telescope (JCMT)^ on Mauna Kea resolves 
the submm EBL into its individual components through 



a combination of unparalleled sensitivity and large field- 
of-vicw. Recent SCUBA surveys (Small, Ivison, & Blain 
1997; Small et al. 1998; Barger et al. 1998; Hughes et al. 
1998; Eales et al. 1999; Lilly et al. 1999; Blain et al. 1999) 
have uncovered a substantial population of dusty galaxies 
with properties similar to those expected for the distant 
counterparts to the most luminous, merging systems ob- 
served locally, the ultraluminous infrared galaxies (ULIGs) 
(Sanders & Mirabel 1996). In this paper we present our 
extensive new wide-area SCUBA surveys of blank regions 
of sky. We analyze our source counts both empirically and 
in terms of semi-analytic models. 

2. OBSERVATIONS 

In Barger et al. (1998) we presented our deepest two 
SCUBA maps of the Lockman Hole and the Hawaii Deep 
Survey field SSA13, the observational details of which can 
be found in that paper. We note, however, that due to 
a bug in the SCUBA observing software during the time 
these data were taken (discovered and fixed 27 July 1998), 
the chopper was started at the updated (Az,El) compo- 
nents corresponding to the requested (RA,Dec) chop throw 
at the beginning of each new 30 min integration but then 
did not update for the remainder of the integration. Con- 
sequently, the off-beams appear smeared in arclets from 
the initial positions. 

Our wide-area SCUBA maps of the Hawaii Deep Survey 
fields SSA17 and SSA22 were taken as mosaics during two 
runs in August 1998 and September 1998 and also during 
an earlier observing shift in October 1997 (totaling 16 ob- 
serving shifts). Our wide-area map of SSA13 was taken 
during a run in February 1999 (7 observing shifts). The 
maps were dithered to prevent any regions of the sky from 
repeatedly falling on bad bolometers. The chop throw dur- 
ing the October 1997 run was azimuthal. During the runs 
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in 1998 and 1999 the chop throw was fixed at a position 
angle of 90 deg so that the negative beams would appear 
45 arcsec on either side East- West of the positive beam. 
The negative images can be restored, thereby increasing 
the effective exposure times on most of the field. Every 
few hours we did a 'skydip' (Lightfoot et al. 1998) to mea- 
sure the zenith atmospheric opacities at 450 and 850-/im, 
and we monitored the 225 GHz sky opacity at all times 
to check for sky stability. Fully sampled wavelength maps 
at 850-/im were completed at median optical depths of 
0.38 during the SSA17 and SSA22 runs and 0.15 during 
the SSA13 run. We performed pointing checks every hour 
during the observations on the blazars 1633-1-382, 2223- 
052, or 1308-1-326. We calibrated our data maps using 
30 arcsec diameter aperture measurements of the positive 
beam in once or twice-nightly beam maps of the primary 
calibration sources Uranus and Mars. 

We reduced all our data in a standard and consistent 
way using the dedicated SCUBA User Reduction Facility 
(SURF; Jenness & Lightfoot 1998). Due to the variation in 
the density of bolometer samples across the maps, there 
is a rapid increase in the noise levels at the very edges; 
thus, we clipped the low exposure edges from our images. 
Our wide-area SCUBA maps and multiwavelength follow- 
up data will be presented elsewhere (Barger & Cowie, in 
preparation). 

3. SOURCE EXTRACTION AND COMPLETENESS 

The SURF reduction routines arbitrarily normalize all 
the data maps in a reduction sequence to the central pixel 
of the first map; thus, the noise levels in a combined im- 
age are determined relative to the quality of the central 
pixel in the first map. To find the absolute noise levels, we 
first eliminated the > 3ct real sources in each field through 



subtraction of an appropriately normalized version of the 
beam profile. We then iteratively adjusted the noise nor- 
malization until the dispersion of the signal to noise ra- 
tio measured at random positions was approximately one. 
This procedure should provide a conservative noise esti- 
mate since it includes both fainter sources and correlated 
noise. 

We scanned the maps at an array of positions separated 
by 6" to locate all > 2.8(7 (to overselect) sources. We note 
that the position determinations in the deep maps are not 
affected by off-beam smearing because it is a symmetrical 
effect. Using the peaked-up positions, we extracted only 
the > 3cr sources. We subtracted the brighter sources from 
the maps before extracting the fainter sources. We cali- 
brated the extracted fluxes, including both positive and 
negative beams, to the 30 arcsec diameter aperture fluxes 
of the brightest sources; the measured fluxes are therefore 
unaffected by the off-beam smearing. 

Our source catalog is given in Table § The noise levels 
for the two deep maps are slightly lower than our initial 
analysis (Barger et al. 1998) because of improvements in 
both the SURF reduction routines and our own extrac- 
tion routines. These improvements result in an increase in 
the number of measured sources for these fields over our 
previously published results. 

To determine the completeness of our source recovery, 
we assumed the shape of the Guiderdoni et al. (1998) 'Sce- 
nario E' model (see §4) for the counts and simulated the ex- 
pected sources by adding appropriately renormalized cali- 
bration sources at random positions to the final maps with 
the real sources removed. We then reran our extraction 
procedure on these images. Any sources detected at the 
> 3(7 level were considered to be recovered. Using this cri- 
terion, the output counts matched the input counts; thus, 
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Fig. 1. — Our 850-/im source counts (solid squares) with Icr error limits (jagged solid lines) are well described by the power- law parame- 
terization in Eq. 1 with a = 0, a = 3.2, and A'o = 3.0 X lO** deg^^ mjy~^ (solid line). The dashed curve shows a smooth extrapolation of our 
fit to EBL measurements using the value a = 0.5. Counts from Blain et al. (1999) (open circles), Hughes et al. (1998) (open triangles), and 
Bales et al. (1999) (open squares) are in good agreement with our data and the empirical fit. Semi-analytic models A and E from Guiderdoni 
et al. (1998) are shown by the dotted curves. 



no incompleteness correction is required. We also investi- 
gated the completeness using the functional form in Eq. 1 
(below) that flattens just below the 2 mJy region and ob- 
tained the same result. At the 3a level, noise may migrate 
sources above and below the threshold, but empirically 
we have found that the net result, given the shape of the 
counts versus flux, is to leave the overall count determina- 
tions unchanged. 

Because of non- uniformity across the maps, the flux sen- 
sitivity levels and area coverages are related. For example, 
at 3cr limiting flux levels of 20, 10, 8, 5, 4, 3, and 2.3 mJy, 
our area coverages, after combining all fields, are 141, 122, 
104, 43, 31, 18, and 7.7 arcmin^, respectively. For com- 
parison, at a 2.1 mJy 4(t flux limit, Hughes et al. (1998) 
have an area coverage of 5.6 arcmin^. The combined data 
of Eales et al. (1999) at 3a flux limits of 2.8, 2.3, and 1.7 
mJy have areas 21.7, 6.7, and 3.4 arcmin^, respectively 
(private communication from S. Eales and S. Lilly); the 
actual la measurement uncertainties given for the 10 and 
14-hour sources in Table 1 of Eales et al. (1999) are mostly 
in the range 0.9-1.1 mJy. 

We present our cumulative 850-/^m source counts per 
square deg (solid squares) with the la errors (jagged solid 
lines) in Fig. 1. The cumulative counts are the sum of the 
inverse areas of all the sources brighter than flux S. We 



can model the differential counts as a function of flux in a 
phenomenological way by fitting 

n(5) =iVo/(a + 5«) (1) 

to the data, where S is measured in mJy. Because the 
data are restricted to the > 2 mJy range, the a param- 
eter makes essentially no difference to a fit of observed 
number versus flux, and we take a = here. A chi- 
square fit gives an optimal index a = 3.2, with a 95 per 
cent confidence range from 2.6 to 3.9, and a normalization 
No = 3.0 X lO'' deg"^ mJy"^ This fit, illustrated in Fig. 1 
with a straight solid line, is consistent with the absence of 
10 mJy sources in our sample. 

Over the measured 2 — 10 mJy range, the contribution 
to the EBL is 9.3 x 10^ mJy deg~^, which corresponds to 
20 to 30 per cent of the total EBL, depending on whether 
we adopt the 850-/im EBL measurement of 4.4 x 10^ mJy 
deg~^ from Fixsen et al. (1998) or 3.1 x 10** mJy deg~^ 
from Puget et al. (1996). Assuming only that the above 
parameterization provides an appropriately smooth con- 
tinuation to fluxes below 2 mJy, we can now fit the EBL 
with a suitable choice of a; for a in the range 0.4 — 1.0 the 
integral of the differential counts over all fluxes matches 
the Fixsen et al. and Puget et al. EBL measurements. 

The shape for a = 0.5 is shown by the dashed line in 
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Fig. 1. The open circles with error bars show the cumu- 
lative counts per square degree from the gravitationally- 
lensed sample of Blain et al. (1999). The Blain et al. sam- 
ple has only one source below 2 mJy (see last two columns 
of their Table 1); at 0.5 mJy we plot their la upper limit. 
The open triangles with error bars show the blank field 
data from Hughes et al. (1998). The open squares show 
the blank field results from Bales et al. (1999); the statis- 
tical uncertainty is only indicated on their brightest flux 
point since the remainder of their data have similar un- 
certainties to those on our data. The agreement between 
all the data sets is excellent and follows the empirical fit 
obtained from our data and the EBL constraint. 

4. DISCUSSION AND CONCLUSIONS 

The IRAS satellite made the first observations of dust 
emission from external galaxies when it measured the sky 
in four bands between 12 and 100-/im. A wide variety of 
infrared luminosities were observed, from normal spirals 
to ULIGs. The spectral energy distributions (SEDs) of 
ULIGs are found to provide reasonable fits to the SEDs 
of submm sources (Ivison et al. 1998). Later in this dis- 
cussion we will assume that the SED of the 'prototypical' 
ULIG Arp 220, the spectral shape of which can be repre- 
sented by a modified blackbody with emissivity and 
a dust temperature of T = 47 K, is appropriate for the 
submm sources. 

It has long been debated whether the dust-enshrouded 
local ULIGs are powered by massive bursts of star for- 
mation induced by violent galaxy-galaxy collisions or by 
AGN activity. A recent mid-infrared spectroscopic survey 
of 15 ULIGs by Genzel et al. (1998) found that 70 - 80 
per cent of the sample are dominantly powered by star for- 
mation and 20 — 30 per cent by a central AGN. Similarly, 
spectroscopic follow-up studies of a gravitationally lensed 
submm sample (Barger et al. 1999; Ivison et al. 1998) in- 
dicate that at least 20 per cent of the sample show some 
AGN activity. We assume in the following discussion that 
a substantial fraction of the submm light arises from star 
formation. 

A variety of models have been proposed to predict the 
submm galaxy population (e.g. Guiderdoni et al. 1998; 
Blain et al. 1999; Eales et al. 1999; Trentham, Blain, & 
Goldader 1999). Guiderdoni et al. (1998) designed a family 
of evolutionary scenarios that ranged from submm source 
populations due solely to the presence of dust in optical 
galaxies (Scenario A) to those due a rapidly increasing 
fraction of ULIGs with redshift (Scenario E). In Fig. 1 we 
compare our blank field 850- /im cumulative source counts 
with these Guiderdoni et al. models (dotted lines). Sce- 
nario A substantially underpredicts the observed 850-/xm 
source counts, whereas Scenario E overpredicts the counts 
by a smaller factor. Because submm fluxes are approxi- 
mately independent of redshift, the predicted number dis- 
tributions versus flux are primarily dependent on the in- 
put luminosity functions. A downward renormalization of 
Scenario E to describe our data would necessitate a steep- 
ening of the luminosity function at the faint end in order 
to match the submm EBL. 

Because of the unique aspects of the submm counts, it 
is also possible to interpret the data in an empirical way. 
With the empirical a — 0.5 parameterization of § ^ we find 
that 60 per cent of the 850-/xm background light arises 



from sources lying between just 0.5 and 2.0 mJy. The 
mean source flux averaged over the whole population is 
0.7 mJy. Because of the steep divergence at low flux levels 
of the empirical fit to the number counts, the turnover of 
the N{> S) distribution will occur in essentially the same 
place independent of the regularization of the divergence 
by specific a-parameter choices; hence, our conclusion that 
the dominant counts occur in the 1 mJy region is not sen- 
sitive to the relative calibration of the EBL and SCUBA 
observations. We estimate the cumulative source density 
in the 1 — 10 mJy range to be 10'' deg~^, which is con- 
sistent with the measurement (0.79 ± 0.3) x 10"* deg~^ by 
Blain et al. (1999); our predicted density for the 0.5 — 2.0 
mJy range is 2.6 x 10^ deg~^. 

Limited and still rather uncertain spectroscopic data 
suggest that most of the submm sources lie in the z = 1 — 3 
redshift range (Barger et al. 1999; Lilly et al. 1999). At 
these redshifts the far-infrared (FIR) luminosity is approx- 
imately independent of the redshift. Thus, if we assume 
an Arp 220-like spectrum with T = 47 K (e.g. Barger et 
al. 1998), the FIR luminosity of a characteristic ^ 1 mJy 
source is in the range 4 — 5 x 10^^ h^^ Lq for a, qo — 0.5 
cosmology (7- 15 x 10" for qo ^ 0.02). The FIR luminos- 
ity provides a measure of the current star formation rate 
(SFR) of massive stars (Scoville & Young 1983; Thronson 
& Telesco 1986), SFR - 1.5 x 10-^° {LfirII^q) Mq yr"!; 
a 1 mJy source would therefore have a star formation rate 
of ^ 70 hgg^ M0 yr~^ for go — 0.5, placing the 'typical' 
submm source at the high end of SFRs in optically-selected 
galaxies (Pettini et al. 1998). If we were to allow the dust 
temperature to go as low as 30 K, Lpm and the corre- 
sponding SFR would be ~ 4 times smaller. 

For a cumulative source density of 4.0 x 10"* deg~^ re- 
quired to reproduce the EBL with 1 mJy sources {{N) = 
EBL/ (5) with (S) - 1 mJy) and redshifts in the 1-3 
range, the average space density is 5 x 10^'^ hgg Mpc"'^ 
for a, qo = 0.5 cosmology (10~^ for qo = 0.02). This space 
density is rather insensitive to the upper cut-off on the 
redshift distribution, dropping by only a factor of ^ 2 or 3 
if we extend the volume calculation to z = 5. For compar- 
ison, the space density of present-day ellipticals is about 
10-3 his Mpc-3 (Marzke et al. 1994; Yoshii & Takahara 
1988). Within the still substantial uncertainty posed by 
the dust temperatures, the estimated star formation rate 
from submm sources is ~ 0.3 hgs Mq yr"^ Mpc"'^ for go = 
0.5, which is nearly an order of magnitude higher than that 
observed in the optical, ^ 0.04 hgs M© yr"^ Mpc"^ (e.g. 
Steidel et al. 1999). 

Several groups (Small et al. 1998; Eales et al. 1999; Lilly 
et al. 1999; Trentham, Blain, & Goldader 1999) have sug- 
gested that the submm sources are associated with major 
merger events giving rise to the formation of spheroidal 
galaxies. The approximate equality of the optical and 
submm backgrounds supports this hypothesis; present-day 
spheroidal and disk populations have roughly comparable 
amounts of metal density, and thus their formation is ex- 
pected to produce comparable amounts of light (Cowie 
1988). 

A star formation rate of 70 hg^ Mq yr"-' Mpc"'^ for 
qo = 0.5 could produce a spheroid of mass 6 x 10^° h^^ Mq 
in 0.8 Gyr. Recent CO observations of two high red- 
shift submm sources by Frayer et al. (1998, 1999) show 
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that each source has enough molecular gas to form an 
L* galaxy. Thus, with a duty cycle that self-consistently 
matches the galaxy number density, the submm popula- 
tion could plausibly evolve into the present-day spheroidal 
population. 

In summary, we have reported new 850-^m source 
counts from wide-area SCUBA imaging data of the Hawaii 
Deep Survey fields SSA13, SSA17, and SSA22 and im- 
proved counts from deep SCUBA imaging data of SSA13 
and the Lockman Hole. Our data rule out models that as- 
sume the submm flux is due only to dust in optical galax- 
ies because such models vastly underpredict the observed 
counts. We find that our cumulative counts, N{> S), 
above a 3cr detection threshold of 2 mJy follow a sim- 
ple power-law behavior. By introducing an additional pa- 



rameter constrained by the EBL, we extrapolate our cu- 
mulative counts below 2 mJy and infer that the bulk of 
the submm EBL resides in sources near 1 mJy, in agree- 
ment with a fluctuation analysis by Hughes et al. (1998) 
and lensed source counts from Blain et al. (1999). The 
submm sources are plausible progenitors of the present- 
day spheroidal population. 

We thank an anonymous referee for helpful suggestions, 
and we thank Stephen Eales and Simon Lilly for commu- 
nicating to us their flux sensitivities and area coverages. 
AJB acknowledges support from NASA through contract 
number P423274 from the University of Arizona, under 
NASA grant NAG5-3042. 
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